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Abstract.  Small  velocity  gradients  in  a  refracting  horizon  have  a  pronounced  effect  on  the 
spectral  amplitudes  of  head  waves.  Negative  velocity  gradients  and  (inelasticity  (Q’1)  result 
in  n  similar  amplitude  decay  with  distance  for  narrow-bandwidth  data.  Positive  velocity 
gradients  result  in  a  net  amplitude  gain  with  distance  compared  with  the  head  wave  from  a 
homogeneous,  perfectly  elastic  refractor.  Wave-theoretical  expressions  for  these  effects  applied 
to  published  amplitude  data  for  the  major  crustal  refraction  branches.  P,  and  P*,  suggest 
that  the  ‘granitic’  crust  in  the  Basin  and  Range  province  has  either  negative  velocity  gradients 
of  the  order  of  10  1  km/sec/km  or  an  anelastic  Q  of  the  order  of  400,  whereas  the  ‘granitic’ 
crust  in  the  eastern  United  States  and  on  the  California  coast  has  slightly  positive  velocity 
gradients.  Similarly,  the  ‘basaltic’  intermediate  layer  appears  to  have  n  negative  gradient  of 
the  order  of  10"*  sec'*  under  the  Snake  River  plain  anil  null  or  slightly  positive  gradients 
under  Lake  Superior  and  Mississippi.  Velocity  gradients  inferred  from  laboratory  measure¬ 
ments  on  granite  and  basic  igneous  rocks,  together  with  published  geothermal  gradients,  are 
generally  consistent  with  the  gradients  inferred  from  amplitude  data. 


1:1  OCITY  GRADIENTS  AND  AMPLITUDES 

cut  wave-theoretical  studies  show  that  the 
spectral  amplitudes  of  critically  refracted  waves 
(head  waves)  are  quite  sensitive  to  small  ve¬ 
locity  gradients  in  the  refracting  horizon  [Cfer- 
ecitfi  anti  Jansky,  1907;  Hill,  1971].  Negative 
gradients  and  anelasticity  (Q'1)  result  in  a  simi¬ 
lar  amplitude  decay  with  distance  for  narrow 
bandwidth  data;  positive  gradients  result  in  a 
net  amplitude  gain  with  distance  with  respect 
to  theoretical  amplitudes  for  head  waves  re¬ 
fracted  from  a  homogeneous  medium  with  in¬ 
finite  Q.  These  theoretical  effects  of  velocity 
gradients  on  head-wave  amplitudes  are  illu¬ 
strated  in  Figure  1.  The  negative  gradient  case 
is  based  on  asymptotic  solutions  obtained  by 
Hill  1 1971  | ;  tin'  positive  gradient  case  is  adapted 
from  ('erectly  and  Jansky  [1967], 

In  a  preliminary  attempt  to  determine  the 
existence  and  distribution  of  velocity  gradients 
in  well-established  crustal  horizons,  these  theo- 

1  Division  of  Geological  and  Planetary  Sciences, 
California  Institute  of  Technology,  Pasadena,  Cali¬ 
fornia,  Contribution  1914.  Publication  authorized 
by  the  Director,  U.  S.  Geological  Survey. 


retical  results  have  been  applied  to  published 
amplitude  data  for  the  major  crustal  travel  time 
branches  P,  and  P*.  In  this  study,  only  first- 
arrival  data  for  these  branches  arc  included. 

RELATION  BETWEEN  O'  AND  VELOCITY  GRADIENTS 

Effective  Q  values  (hereafter  referred  to  as 
Q')  were  computed  for  each  amplitude  data  set 
by  determining  in  a  least  squares  sense  the  net 
gain  (  —  Q')  or  decay  (  +  Q')  of  the  observed 
data  with  respect  to  theoretical  head-wave 
amplitudes  from  homogeneous,  infinite  Q  refrac¬ 
tors.  The  results  are  interpreted  in  terms  of 
negative  velocity  gradients  for  Q'  values  smaller 
than  10*,  null  gradients  for  Q'  values  greater 
than  10s,  and  positive  gradients  for  negative  Q' 
values  :iu  summarized  in  Table  1.  Four  examples 
of  the  amplitude  data  showing  the  least  squares 
fit  for  Q'  are  given  in  Figure  2. 

O'  VERSUS  REGIONAL  HEAT  FLOW 

The  distribution  of  Q'  values  for  the  data  con¬ 
sidered  is  shown  in  Figure  3,  together  with  re¬ 
gional  heat  flow  contours  presented  by  Archam- 
beau  et  al.  [1968],  From  the  figure,  we  see  that 
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Fig.  1.  Amplitude  curves  illustrating  cfTect  of 
positive  ( /3  =  10'*  sec"’),  null  (ft  —  t>) .  mid  nega- 
tivc  (ft  =  — 10‘*  see"1)  velocity  gradients  on  lie.nl- 
wnvc  amplitudes,  together  with  effective  (}  values 
(Q')  with  respect  to  the  ft  —  0,  Q  =  oc  case.  The 
reflected  wave  amplitude  r.  which  is  insmsitive  to 
small  gradients  in  the  refracting  medium,  is  shown 
for  reference.  The  arrow  at  C.I).  indicates  the  ray- 
theoretical  critical  point.  These  curves  were  gen¬ 
erated  for  G-IIz  waves  in  a  (1.4-km/see  layer  .'10 
km  thick  over  an  8.0-km/sec  half-space.  The  posi¬ 
tive  gradient  ease  is  adapted  from  Ccrvcity  ami 
Jansky  [19671. 

there  is  a  general  inverse  correlation  between 
Q'  values  and  heat  flow.  Q'  values  in  the  'normal’ 
heat  flow  provinces  of  the  eastern  United  States 
and  west  coast  arc  greater  than  1(1*  or  are  nega¬ 
tive;  these  values  imply  null  or  small  positive 
velocity  gradients  for  both  the  upper  portions 
of  the  crystalline  crust  (‘granitic  layer’)  and  the 
intermediate  layer  in  these  regions.  Q'  values  in 
the  Basin  and  Range  high  heat-flow  province, 
which  probably  includes  the  Snake  River  plain, 
are  generally  smaller  than  10*.  These  low  Q' 
values  common  to  the  Basin  and  Range  province 
may  be  due  to:  1,  scattering  of  the  critically 
refracted  waves  by  relief  on  the  refractor  as¬ 
sociated  with  Basin  and  Range  faulting;  2,  a 
temperature-dependent  .inelastic  Q ;  or  3,  nega¬ 
tive  velocity  gradients  in  the  upper  portions  of 
the  crystalline  crust  associated  with  high  geo¬ 
thermal  gradients. 

Available  information  is  insufficient  to  quan¬ 
titatively  assess  effects  of  the  first  two  factors, 
and  they  remain  as  possible  contributing  factors. 
The  third  factor  can  be  assessed  quantitatively 
and  is  found  to  be  consistent  with  the  correla¬ 
tion  between  Q'  and  heat  flow  mentioned  above, 
lleat  flow  studies  suggest  that  the  crustal  geo¬ 


thermal  gradients  in  the  Basin  and  Range  heat 
flow  province  is  ,'it)°/kin  or  possibly  somewhat 
higher  [  Hay  et  al ,  10(iS;  Laehenhrueh,  1070; 
Minster  (iinl  Arehmnbemi ,  1*170 J .  This  geother¬ 
mal  gradient,  together  with  partial  derivative 
data  for  /'-wave  velocities  in  granites  with  re¬ 
spect  to  pressure,  (ill’.,  i)/’)r,  and  temperature, 
(dr.,  (V/’l,..  measured  in  the  laboratory  \  Ilnghcs 
(iml  Mniireth,  I!(5t>|.  suggests  negative  velocity 
gradients  of 

d  ~  -0.S  X  10  2  km/sec/km 

in  tin'  ti.0-km  sec  ‘granitic’  horizon  of  the  Basin 
and  Range  province.  'Pile  lower  geothermal  gra¬ 
dients  (about  I.iVkm)  that  arc  associated  with 
the  ‘normal’  heat  flow  province,  when  combined 
with  tin*  same  partial  derivative  data,  suggest 
positive  velocity  gradients  of 

d  ~  0.5  X  10  2  km/sec/km 

in  the  upper  crystalline  crust.  Thus,  the  differ¬ 
ences  in  velocity  gradients  in  the  upper  part  of 
the  crystalline  crust  (/',  refractor)  between  the 
Basin  and  Range  and  eastern  United  Sta!es- 
west  coast  heat  flow  provinces  inferred  sepa¬ 
rately  from  /'„  amplitudes  and  geothermal  gradi¬ 
ents  are  reasonably  consistent. 

Similar  calculations  using  the  partial  deriva¬ 
tive  data  for  /'  waves  in  basic  igneous  rocks  re¬ 
ported  by  Hughes  and  Maurette  [1057]  indicate 
that  the  ‘basaltic’  ti.7-km/sec  intermediate  layer 
has  negative  velocity  gradients  of 

d  ~  —  1.6X  10  '' km/sec/km 
beneath  the  Snake  River  plain  and 

d  ~  —0  3  X  10  "  km/sec/km 

beneath  Lake  Superior  and  Mississippi.  As  be¬ 
fore,  a  geothermal  gradient  of  30°/l;m  is  as¬ 
sumed  for  the  high  heat  flow  region  (the  Snake 
River  plain)  and  15°/km  is  assumed  for  the 
‘normal’  heat  flow  regions  (Lake  Superior  and 
Mississippi).  These  velocity  gradients  inferred 
from  geothermal  data  are  more  negative  than 
th°se  implied  by  the  Q'  values  obtained  from 
/'*  amplitude  data  in  these  regions  (Table  1). 
If  we  accept  this  result  at  face  value,  then  a 
compositional  gradient  (increasingly  mafic  with 
depth)  is  required  in  the  upper  parts  of  these 
intermediate  layers  to  bring  the  velocity  gradi¬ 
ent  estimated  from  the  geothermal  gradient 
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TABU'!  1.  Summary  of  Q'  and  Velocity  (iradients  from  Crustal  Amplitudes  Data 
Standard  error  refers  to  Q'~x  fit  to  amplit  tide  data,  d  is  velocity  gradient,  and  'null'  and  d  indicate 
zero  and  positive  gradients,  respectively. 


Region  and  Profile 

Phase 

tf 

(J 

Std.  Error 
in  Q' 

d(sec  *)" 

Source 

Basin  and  ltange 
Fallon- Eureka 

1\ 

■m 

2.10 

0.91 

-15 

Eaton  [1963] 

Fnllon-.S.F. 

p. 

760 

1 .22 

0.36 

-  9.7 

Eaton  (i!)63| 

Fallon-Owens  V. 

p. 

-446 

—  2  2 

1.16 

+ 

Eaton  [1963] 

Eurcka-Fnllon 

p. 

972 

I  .03 

0.47 

-  7.5 

Eaton  (1963| 

Eureka-North 

- 1290 

-0.77 

0.88 

+ 

Hill  A*  Pukixer  [19061 

Ml.  City-South 

403 

2.18 

3.71 

-18 

Hill  it  Pakixer  (I960) 

NTS- East 

p. 

117 

8.50 

0.87 

-26 

Hyall  it  Stuart  (19631 

California 

S.F.-Fallou 

3810 

0 . 262 

0.076 

nidi 

Eaton  (1963) 

S.F.-S.  Monica 

p. 

- 1580 

-0.03 

1.13 

-f 

Healy  (1963) 

Camp  Roberts 

i\ 

2860 

0.349 

0.87 

null 

Hcaly  (1963) 

S.  Monica-L.  Mead 

i\ 

230 

4.34 

2.17 

-15 

Roller  it  Healy  (1963) 

San  Juan  (0.06 ) 

p. 

-47 

-21. 

4.9 

-F 

Stewart  (19685) 

San  Juan  (6.35) 

/’? 

54 

18. 

3.7 

-20 

Stewart  19685 

Colorado  Plateau 

Hanksvillc 

1260 

0.795 

1 .02 

null 

Roller  (1965) 

Chiide 

p. 

221 

4.53 

1.33 

-16 

Roller  -1965] 

Missouri 

Hannibal 

r. 

-613 

-1.63 

0.86 

+ 

Stewart  (1968a) 

Swan  L.-lIannibal 

p. 

-549 

-1.82 

0.54 

+ 

Stewart  [1968a| 

Swan  L.-St.  Joseph 

p. 

- 1960 

-0.511 

0.633 

+ 

Stewart  (1968a) 

St.  Joseph 

p ; 

-475 

-2.10 

0.63 

+ 

Stewart  (1968a) 

Mississippi 

/*• 

-382 

-2.62 

0.92 

+ 

Warren  el  al.  1 1966) 

Lake  Superior 

r* 

31,700 

0.0315 

0.604 

null 

O'Brien  (1968) 

Snake  River  plain 

Boise-South 

r* 

337 

2.97 

0.84 

-  4.4 

Hill  it  Pakiser  (1966) 

•  1  X  io-‘. 
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Fig,  2.  Examples  of  P,  amplitude  data  in  four  regions:  a,  Fallon  toward  Eureka  in  the 
Basin  and  Range  [Eaton,  1963];  b,  San  Francisco  cast  in  California  [Eaton,  1963];  c,  Hanks- 
villo  north  in  (he  Colorado  1  lateau  [Roller,  1965] ;  and  d,  Hannibal  west  in  Missouri  [Stewart, 
1968fiJ.  Continuous  line  is  least  squares  fit  for  Q'  through  the  amplitude  data.  Broken  line 
shows  amplitude  decay  with  Q’  =  oo  for  comparison. 


i 

i 

j 


i 

I 


71 


CONTINENTAL  SElSMULUtn 


Fig.  3.  Locatin':  of  seismic  profiles  and  the  relation  of  Q'  to  heat  flow.  Number  beside  each 
profile  is  Q'  X  10  3.  The  dark  broken  contour  line  separates  the  region  in  which  Q’  is  less 
than  10”  from  that  in  which  it  is  either  greater  than  10*  or  negative.  Regional  heat  flow  pat¬ 
terns  are  adapted  from  Archnmbcnu  et  al.  '19681.  Numbers  associated  with  heat  flow  contours 
have  units  Mcal/cnr'/sec. 


(eg.,  —It»  X  10’  sec'1  in  the  Snake  River 
plain)  to  the  level  of  the  velocity  gradient  esti¬ 
mated  from  ()’  (4  X  10  *  sec'1  in  the  Snake 
River  plain).  However,  neither  the  amplitude 
nor  the  thermal  data  arc  of  sufficient  accuracy 
to  provide  much  confidence  for  this  suggested 
compositional  gradient  in  the  intermediate  lay¬ 
ers. 

CONCLUSIONS 

Thus,  on  the  basis  of  P,  amplitudes  and  ther¬ 
mal  data,  \vc  conclude  that  a  slight  P-wave  low 
velocity  zone  may  exist  in  the  upper  crystalline 
crust  in  the  Rasin  and  Range  high  heat  flow 
province  and  that  such  zones  arc  unlikely  in  the 
eastern  United  States  and  west  coast  normal 
heat  flow  provinces.  In  the  Basin  and  Range 
province,  the  crustal  low  velocity  zone  would 
have  the  form  of  a  gradual  decrease  in  velocity 
from  the  top  of  the  crystalline  crust  downward 
(at  a  maximum  rate  of  about  1  X  10'*  km/ 
sec  km)  and  would  terminate  rather  abruptly 
at  the  top  of  the  intermediate  layer.  The  abrupt, 
pronounced  crustal  low  velocity  zone  at  depths 
of  about  It)  km  of  the  type  proposed  by  Mueller 


and  Landisman  flfttiOj  could  be  present  in 
cither  the  Basin  and  Range  province  or  the 
eastern  United  States,  but  such  a  zone  would  be 
difficult  to  detect  using  the  data  and  methods 
described  in  this  paper. 
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DISCUSSION 


Porath:  You  have  a  uniform  f).4-km/sec  ve¬ 
locity  for  the  Columbia  plateau.  Might  you  not 
expect  a  lower  velocity  at  depth  because  the 
plateau  basalts  probably  cover  an  ancient  gra¬ 
nitic  crust  ? 

Hill :  The  data  only  give  the  time  it  takes  a 
Pvmv  in  ir-iVel  ft  in  Hr  Yt  j!ij  hi  llm  svtrhor’ 
there  is  no  way  of  dividing  the  crust  up  into 
layers,  so  this  is  just  an  average  velocity. 

Mei  er  •  J  )W  VHLI  iUlfchLT  UUItg  h.  fjlWrf  t) 
and  positive  gradient  in  the  Lake  Superior  re¬ 
gion  0 


Hill:  All  I  did  was  to  fit  the  homogeneous 
head-wave  potential,  derived  from  the  ampli¬ 
tude  data  published  by  O’Brien,  and  make  a 
least  square  solution  for  Q.  This  gives  an  esti¬ 
mate  of  the  lower  limit  for  the  gradient  and  an 
upper  limit  for  Q. 

•  What  ifl  the  pTijvirnl  nignlffrnnre  of 
a  negative  Q ? 

Hill:  That  is  an  artifact;  in  reality  there 
tun  ill  l*1  i  [iL/i'iir“  Q  sru!  „  pnsitiv-  velocity  gra¬ 
dient  to  give  the  amplitude  increase. 
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